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ABSTRACT 

We have been able to connect the statistics of the observed double image gravitational 
lenses to the general properties of the internal structure of dark matter haloes. Our 
analytical theory for the GNFW lenses with parametrized cusp slope (a) gives us a 
relation connecting the cusp slope of the lcnsing profile to the observed magnification 
ratio of the produced images and location of the optical axis. The relation does not 
depend on cosmology, total lens mass, concentration or redshifts of the the lens and 
the lensed object. Simple geometry of axially symmetric lensing and aforementioned 
relation enables us to define a threshold value a CSL for the cusp slope, independent 
from location of the optical axis. The threshold cusp slope value a = a C sL is the 
shallowest slope for the inner part of the GNFW profile that can produce the ob- 
served magnification ratio with any lensing configuration. We use distribution of these 
threshold values in a statistical study of the double image lenses in order to limit the 
possible cusp slope values, and identify whether there exists a population of haloes 
with similar profiles. Our theoretical fit indicates that within our sample of double 
image gravitational lenses, most of the haloes have cusp slope a = —1.95 ± 0.02. We 
have also found an indication of a second population of lenses with a cusp slope value 
a = —1.49 ± 0.09. We estimate that there is about 99 per cent probability that the 
observed feature in the threshold value limit distribution is produced by the second 
population of lenses, with their own characteristic density profile. The data indicating 
the exact characteristics of the sub-population is noisy. Roughly one out of six haloes 
within the sample belong to this shallower cusp slope group. We investigate errors in 
our analysis by constructing mock catalogues with the Monte-Carlo method. 

Key words: gravitational lensing, relativity, methods: analytical, galaxies: haloes, 
cosmology: dark matter 



1 INTRODUCTION 

The large scale structure formation in the universe is 
dominated by the gravitational evolution o f dark matter 
(DM) and the expansio n of the universe (|Peeblesl 1 19821 : 
iBlumenthal et al. Ifl983 : iDavis et al1ll985l ). Currently, the 
best candidate for the DM component is a weakly in- 
teracting massive particle (WIMP). A host of candidate 
particles can be suggested from the various extensions of 
the standard particle theory, fo r example supersymmetry 
Ellis et al.lll99ll; iMartin 1119981) or Kalu za-Klein concepts 



Cheng et al.ll2002l : IServant fc Taitll2003l ). 



Theories concerning structure formation usually pos- 
tulate non-interacting dark matter that is affected only 
by gravity. Since the last decade, a large number of 
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studies employing analytical theories and/or N-body sim- 
ulations have been conducted on this basis. There are 
also several variants of this approach, in which DM 
has been attributed with additional properties, such as 
being s e lf-inte r acting (e.g. Spergel fc Steindhardt I 2000l; 
Burkertj |200C|; iRomeel et all l200ll: iD'Onghia fc Burkertl 



20031 : ID'Onghia, Firmani fc Chincarinil |2003|. and refer- 



ences therein), decaying (e.g. 
iBorzumati. Bringmann fc Ulliol 

lated equation of state (for example Austin et al. 120081 ) 



Adbelqader fc Melial 120081 : 
2008|) or having a po stu- 



At the linear regime, the current state-of-the-art cos- 
mological theory, ACDM-model with inflation, is very suc- 
cessful and can explain most of the observations. How- 
ever, problems arise at the nonlinear regime, where the 
DM halo is decoupled from the expansion of the universe 
and evolves through self-gravitation and gravitational in- 
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teractions with it s environment (|Moore et al.l [l998. 1999; 
iPower et al1l2003l ). 

At the high density region of dark matter halo cores, 
the physical properties of the dark matter particles and 
their interactions shoul d beco m e important. A ccording to 
ISpereel fc Steindhardt I l|2000l ): iBurkertl £2000) the weak- 
interaction cross section of the dark matter particle can have 
significant effect on the shape of the density profile on cos- 
mological ti me scales. The c ross s ection has also been con- 
strained by iRandall et al~l (I2002TI using the Bullet galaxy 
cluster. On the other hand. iBeacom. Bell fc Mack I (|2007n 
determine upper limit for the weak-interaction cross section 
from the cosmic diffuse neutrino background, and argue that 
dark matter haloes cannot be significantly altered by dark 
matter particle annihilations. It is clear, that observations 
and theoretical studies of dark matter haloes can help us 
narrow down the properties of the candidate particles. 

The study of dark matter halo format ion at the non- 
linear regime by iNavarro et alJ (|l996l . 119971 ) found that the 
dark matter haloes follow roughly a universal radially sym- 
metric density profile (NFW). In the study, one halo con- 
sisted of several thousand gravitating DM-particles. Sim- 
ilar N-body experiments have been conducted regularly 
with increasing number of DM- particles and better resolu - 
tion by sever a l auth ors, see e.g. iMoore et al ] l|l99Sl . ll999h : 
iGhigna et ail (|200Cj ). 

Since then, several N-body study inspired models for 
the DM-haloes have been proposed. They all share some 
common properties; at the outer fringes, the halo follows 
p oc r~ 3 profile. After the transition region at scale radius 
r s , the profile is changed to value p oc r a representing a 
cusped (or flat) core. 

An ad hoc explanation for this behaviour is that the 
outer region is in the state of 'inflow', where the dark mat- 
ter decouples from the general expansion of the universe, and 
streams towards the core region of the halo. This is charac- 
terized by dominating radial component of the velocity field. 
Inside the scale radius r B , where the velocity field is more 
thermalized, the halo is composed of captured dark mat- 
ter component that has passed thro ugh the core and turned 
back a t l east once. See for example iDehnen fc McLaughlin! 
|2005l ) or lHansenl (120091) for further details. 

In their study. INavarro et all l| 19961 , 1 19971 ) found an uni- 
versal cusp slope a = — 1, valid for haloes at extensive range 
of size scales. More recent N-body studies have resulted an 
array of different slope values ranging from a — —1 to 
a = — 2. This ambiguity in resolving the cuspiness of the 
DM haloes with N-body simulations, and contradicting ob- 
servations from galactic dynamics has been one of the central 
weaknesses of the ACDM paradigm. 

The observed shape of the density profile can be mea- 
sured indirectly by analyzing internal dynamics and kine- 
matics of nearby galaxies, measuring X-ray temperature pro- 
file, using Sunyaev-Zel'dovich effect, or with weak lensing 
measurements of galaxy clusters. Several authors has pub- 
lished results characterizing the pro file shape by combining 
some of these observations, see e.g. iMahdavi et"ai~l (|2007l ) 
and references therein. 

These studies either assume a profile shape in 
advance or try to acquire a fit for it. In both cases, 
the N FW pro fil es (a ~ —1.0) have proven success- 
ful (G avazzil 120051 ; IPointecouteau. Arnaud fc Prat j 



120051 ; IVikhlinin et~aH 120061 : iLimousin et"aH 120081 ). al- 
though shallower density profiles a ~ —0. 5 have 
been reported |Sanderson. Finoguenov fc Mohrl 120051 ; 



IVoigt fc Fa bian 2006 ). On the other hand, lensing stud- 
ies by Sereno. Lubini fc Jetzerl (j2009) estimated slightly 
steeper density profile for the cluster AC114 than the 
canonical NFW value. 

Rotation curves and analysis of dynamics for late type 
spiral galaxies indicate that they must have shallower n on- 
cusped cores (jValenzuela et al ] |2007l : lDonato et alj[2009h . If 
a galaxy has a bar, it will interact with the halo and thi s 
should slow-down the bar rotation, see e.g. IWeinberdl|l985h . 
It has been claimed that bars rotate so fast (corotation reso- 
nance close to bar radius) that they could not be embedded 
into strongly concentrated halos. On the other hand, the in- 
teraction between the bar an d halo have been suggested t o 
make the cusp more shallow (jDebattista fc Sellw ood 2000). 
However, there seems to be possible issues with the reso- 
lution of N-bod y models l|Weinberg fc K atz 2007]), see also 
ISellwoodl i|200Sf ). 

Interestingly, in recent high resolution N-body models 
bv lDubinski. Berentzen fc Shlosmanl (|2009l ) . the bar stayed 
fast and did not destroy the cusp. Further complications 
are provided by rec ent results indicating that some small 
bars rotate slowly llRautiainen. Salo fc Laurikainenl 120081 : 



Ch cmm fc Hernandez! 120091). Slow rotation for such bars 
was suggested already bv lElmegreen fc Elmegreenl (|l985l ) on 
morphological grounds, so it is possible that their origin is 
somehow different from larger bars that are studied more 
frequently. 

The dwarf galaxies at low surface brightness are consid- 
ered as being dominated by dark matter component. How- 
ever, direct measurements of the rotation curves of these 
objects indicat e a constant density core w ithout a cusp, see 
IBurkertl (|l995l ) or lZackrisson et al. I |2006h . and correspond- 
ing references therein. 

The newest development in the constraints for the dark 
matter halo profiles comes from the observational parti- 
cle physics. The WMAP probe has observed a microwave 
"haze" around the Galactic core. This can be explained by 
hypothetical dark matter particles annihilating at the high 
density region of the Milky Way core. The "haze" is inter- 
preted as emitted synchrotron radiation from the produced 
electron-positron pairs, and it matches predicted emission 
from a dark matter halo wit h a cusped profile, p oc r~ 12 
withi n the inner kiloparsecs l|Hooper. Finkbeiner fc Doblerl 
2007). Since there is no definite direct detection of the dark 
matter particle to date, this can be considered as a hypo- 
thetical scenario yet to be proved or disproved. 

The largest coherent sample of the galaxy scale strong 
gravitational lenses has b een measured with t he HST in the 
Sloan Lens ACS survey (|Bolton et al.ll2006h .l Gavazzi et al.l 
(2007) made a statistical study about subset of these lenses, 
and confirmed that on average, an isothermal profile a ~ 
—2 fits the lensing profiles. Several strong gravitational 
lensing studies have also suggested shallo wer cusp slopes 
l|Sand. Treu fc Ellisll2002l ; ISand et al.ll2004l ) at a ~ -0.5. 

However, flat profiles hav e severe problems in stron g 
lensing. In our previous paper [Mutka fc Mahonenl (2006), 
we developed a semi-analytical theory for axially symmetric 
gravitational lensing and emplo yed it to the generalizatio n 
of the NFW profile (GNFW), bv lNavarro et all l| 19961 . Il997h ; 
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Strong lensing cusp slope limit 




inner caustic image should be very common. The estimated 
probability for a triple image lens configuration is 



0.5 



Figure 1. The strong lensing condition, equation J2} is function 
of the mass M of the lens. This plot shows the required mass 
for the lens as a function of the strong lensing cusp slope limit 
c»sl. The strong lensing is not possible on the right hand side 
(below) of each curve. The curves with concentration parameter 
values ci = 5, 10, 15, 20, 25, 30 are plotted for the lens at redshift 
z = 0.3, and the source at redshift 2 = 1.0. The dashed line is a 
numerical version of the condition approximated by equation J2) 
as prescribed in the appendix lAl 



IZhaol (|l996h . In the GNFW profile, value of the cusp slope 
is a free parameter: 



p oc r a (l + r)~ 



(1) 



Note that we use convention in which the cusp slope a has 
negative values. In axial lensing context, the GNFW halos 
with cusp slope sh allower than a ~ — 1 have severe problems 
in strong lensing (|Mutka fc Mahonenll200fJ ). 

The cusp slope limit for strong lensing a S]j is a function 
of the cosmological parameters, lensing geometry, lens mass 
and lens concentration. Basically, at this limit the radius of 
the Einstein ring for a GNFW lens goes to zero, after which 
the strong lensing configurations become impossible. A good 
approximation for this limit is 



3 1 L 
a SL = -- + -Jl 



1600r 3 p cr cf(l + ci) 
3E cr [(l + ci)log(l + ci) + Cl ] 



(2) 



that produces values slightly above the a ~ — 1, with reason- 
able lens masses and concentrations. Figure [T] shows the re- 
quired lens mass as function of the strong lensing condition, 
calculated with equation ((2} and corresponding fully numer- 
ical solution. With usual lens masses order of M ~ 10 12 M Q , 
the strong lensing is limited close to a ~ — 1. Only excessive 
concentration allows lens systems with slope a ~ —0.5 that 
can produce multiple images. Appendix [X] presents defini- 
tions of r„, pa and £„ and derivation for the equation ([2]). 

It can also be shown that, if dark matter haloes have 
shallower profile (a —1.5), triple image lenses with visible 



Ptriple(a) 



a 

T 



(3) 



and there are only three observed lens systems with (sus- 
pected) visible inner caustic image out of one hundred. This 
gives an average cusp slope value a ~ —1.97, which is closer 
to the isothermal value. See section [2~T1 for derivation of the 
equation (|3]) and further details. Note that equations ((2)1 and 
((3} are also valid with reasonable values of eccentricity for 
elliptic lens haloes. See section [2] for further discussion. 

We investigate the characteristics of the lensing profiles 
by deriving a cusp slope limit (CSL) value Ocsl for the shal- 
lowest cusp slope, that can produce the required flux differ- 
ence for the images with any lensing geometry (i.e. it must 
hold a < «csl for the lensing profile). We assume random 
alignments for the sources, and use the CSL value in sta- 
tistical analysis of the double image lenses. By examining 
distribution of these threshold values, we can characterize 
the general properties of lensing profiles. For this purpose, 
we composed a catalogue of the double image lenses exhibit- 
ing properties of radial symmetry. Our sample is based on 
the known cases of the lensed quasars listed in CASTLES 
survey. 

Our method relies on the observed flux ratios of double 
image lenses. Because we use image flux ratios, information 
on the lens and source redshifts, lens mass and concentra- 
tion, and cosmology is not required. With our formulation, 
the intrinsic halo properties are separated from these quanti- 
ties. Therefore, the statistical method does not require the- 
oretical distributions characterizing lens and source popu- 
lations (such as Schechter function coupled with the mass- 
luminosity relation for lenses, or quasar luminosity functions 
for sources), which makes it possible to avoid number of dif- 
ferent uncertainties in the analysis. 

However, there are several factors that can corrupt the 
observed fluxes from the predicted magnification produced 
by the lens. We examine effects by the lens ellipticities, the 
substructure in the lensing halo and time-delay from a vari- 
able source. This is done by constructing mock lens cata- 
logues with the Monte Carlo method when accounting for 
these factors. 

The next section presents a brief review of our lens- 
ing theory, formulation for the CSL value, triple image lens 
probability and description for theoretical CSL distribution. 
The third section describes the selected lens systems and the 
catalogue used in the study. The Monte Carlo testing of the 
CSL statistics is presented in section four. Finally, we em- 
ploy the statistical CSL analysis on real double image lens 
data in section five. The last section contains conclusions 
and discussion regarding the study. The derivation for the 
cusp slope limit for strong lensing is presented in appendix 

in 



2 THEORY 

Here we briefly revie w the relevant lensing theor y presented 
in our previous paper l|Mutka fc MahSn en 120061 ) , and derive 
the cusp slope limit (CSL) value. The CSL value is based on 
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lens plane 



optical axis 




Figure 2. The lens equation normalized with the Einstein radius 
(J4j produces three images fei, &2 and k^. The first image k\ is 
located between the Einstein radius (Z/q = 0orfc = l) and the 
first critical radius k = v. The second image is at the opposite side 
of the optical axis, between the Einstein radius and the second 
critical radius, equation J5} (the largest possible distance from 
the optical axis for the multiple imaged source). The third image 
fc3 (inner image) is inside the first critical radius, and it is usually 
strongly demagnified. With large values of k, outside the second 
critical radius, the lensing curve should approach to l/q = k. The 
chosen lensing curve strongly exaggerates size of the triple image 
region. 



the analytical formulation of the magnification ratio of the 
lensed images, and its solutions. 

In gravitational lensing, the source image is lensed by 
the gravitational field of the lens object. The source image 
that is thought to be confined on a plane at the source dis- 
tance D s , produces observed images on the lens plane (at 
the lens object distance D L ). Mapping from the source plane 
to the lens plane is described by the lens equation. 

In axially symmetric lensing, we choose the line con- 
necting the observer and the lens as origin for the frame 
of reference. Corresponding radial coordinate at the source 
plane is denoted with y and at t he lens plane with x, as in 
ISchneider. Enters fc Falcol l| 19921 ). 

When the source is located at the origin (y = 0), it is 
mapped to the lens plane as an Einstein ring, that has radius 
x = x e . The axially symmetric lens equation for the GNFW 
density profile, normalized with the Einstein radius x e , can 
be written as 



qk(l- 

fklog 



\k\ a+1 ) 
\k\ 



if a / -1, -2 < a <0 
if a = -1 



(4) 



where the radial coordinate at the source plane is I = 
y/(x B r s ) and corresponding radial coordinate at the lens 
plane is k — x/x c = r/(r B x Q ). Scale radius r s , and con- 
stants q and /io depend on the cosmology, the mass and the 
concentration of the lens object, and the geometrical setup 
of the lensing. Note that our cusp slope value is a negative 
number. The exact definition for the scale radius r s and the 
corresponding unnorm alized lens equation c an be found in 
appendix 1X1 and paper iMutka fc Mahonen I (|200fJ ). 

With this notation, the Einstein radius is located at 
k = 1. The first critical radius at fc crl = v = (a + 2) _1/(a+1) 
and the corresponding twin image at the opposite side of 
the lens has a good approximation 



k CT 2 — 



i + yi + 2^ 

a + 2 



(5) 



The critical radii (k — 1 and k — fc orl = v) are connected 
to the caust ic curves where the magnifi cation diverges, see 
figure H or IMutka fc Mahonen I |2006) for further details. 
Note that although the image at k cll is at the critical curve, 
the corresponding twin image at fc cr2 is not. 

The approximation in the lens equation @ breaks down 
with large values of the radial coordinate k, and therefore 
the lens equation Q does not behave asymptotically as it 
should. This asymptotic behaviour can be corrected by writ- 
ing the lens equation (a/-l) piece- wisely as 



I = 



when k $J kt 
when k > k B 



gfc(i-ifcr +l ), 

where the lens equation (0 is divided at 
'2(q ^\VC«+l) 



11 

q(a + 3) 



(6) 



(7) 



in order to avoid negative surface densities. 

A lensed image at location k ^ fc B is magnified by the 
lens as 



/'• : 



-2 I 1 

q i 



|fc| QH 



{a + 2)\k 



ct+l I - 1 



(8) 



Note, that the factor q dependency is contained into a co- 
efficient multiplying the magnification. Therefore, the q de- 
pendency is cancelled out when considering ratio of image 
magnifications in strong lensing cases. 

The positive definitive lensed images at ki, ki and ks 
that are solutions of the lens equation Q or (JSJl , are labelled 
as in figure [2] The image closest to the source, that has 
positive parity, is located at fe. The images at fci and k$ are 
at the opposite side of the lens with negative and positive 
parities (the image at k$ is usually strongly demagnified). 
The solutions ki and k^ of a source at I in the lens equation 
must satisfy 



a + 2 _|_ k a+2 



ki + k2 



— 1, when a/-l. 



(9) 



Now we can parameterize the lens equation with coordinate 
ratio 8 = k\jki = x\/x2 = 61/62 ^ 1, where 6\ and 62 are 
observed angular distances of the images from the optical 
axis on the sky correspondingly. The advantage of the pa- 
rameter 8 is that it is directly observable, if the location of 
the optical axis is known. 

The ratio of the image coordinates 6 = fci/fo has a 
minimum value at 
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v{a + 2) 



:+ 1 + 



(10) 



CSL solution for M = 0.7 



corresponding to the most asymmetric configuration for the 
strong lensing for a source image at / = Z max . If the image 
coordinate ratio 9 is decreased below this value, it describes 
the ratio 6 = k$/k2 for the demagnified inner caustic image. 
Thus, # min is also the smallest value for the 9 — ki/ki- 

When relation (0 is employed, and -1, the image 
coordinates become 



ki 
fca 



8 [lT9^] 

r 1+9. 1" 



!/(«+!) 



(11) 



and the lens equation can be written as 

l/Cc+l) 



l + l 



1 + 



Qa + 2 



1 



k< k B . (12) 



a + 2 _|_ I ' 

Now the magnification ratio M of the lensed images 
as a function of the cusp slope a and ratio of the image 
coordinates 6 = fci/fe becomes 

EL 

fl2 



M{a,9) 



|0 Q+2 + l-0 Q! + 1 (l + 0)(a + 2)|' r ' v ; 

as presented in iMutka fc Mahonenl (120061 ). Note that the 
derived relation (|13p between the location of the optical axis, 
the image magnification ratio and the value of the cusp slope 
does not depend on the mass or the concentration of the 
lens, the cosmological model, or the redshifts of the source 
and lens object. This makes it an ideal tool for studying the 
structure of the lensing halos by separating their intrinsic 
properties from the extrinsic lensing conditions. 

When the magnification (flux) ratio of the observed im- 
ages (Mo) is measured, the location of the optical axis (0) 
as a function of the cusp slope parameter a can be solved 
from equation (|13|l . Basically, there are two different groups 
of solutions corresponding to triple and double image lens 
configurations. 

If the solutions are plotted at (a, 6*)-plane, the double 
image solutions vanish at a — lycsl and 8 = #csl, where 

» = "CSL 
a = a CSL 



This means that the double image solutions are restricted 
to cusp slopes a G [—2, «csl]- But on the other hand, using 
the relation (|13|) . this point can be related to the observed 
magnification ratio with 



M(q csl ,0 C sl) = M . 



(15) 

Pair of equations (I14|l and (|15[) . can be used to solve #csl 
and «csl for a lens system with observed magnification ratio 
Mo — Hi/fi2 = F1/F2 (where F\ and F2 are corresponding 
observed fluxes for the images) . See figure [3] 

At the configuration (qcsl,#csl), the lens has just 
enough 'power' to produce the magnification difference for 
the observed flux ratio. For shallower cusp slope values 
a > cvcsl the observed flux ratio cannot be reached with 
any double image configuration. 

The absolute maximum for the possible values of the 
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Figure 3. The CSL solution for a double image lens system 
with image flux ratio M(a, 8) = Mo = 0.7, as calculated from 
equations 1 ] 1 4[ ) and 115> . Horizontal axis shows coordinate ratio 
9 = ki/ki and vertical axis profile slope a. Curve dM/d£> = 
intersects the magnification ratio solutions M(a, 0) = Mo at the 
shallowest possible cusp slope value a = ocsl- The shaded area 
corresponds to the inner image solutions, where k\ is replaced by 
&3, and triple image solutions are limited nearby 6 = m i„. 



cusp slope a can be calculated from the measured flux ra- 
tio for each observed double image lens system. This is the 
previously discussed cusp slope limit, abbreviated as CSL. 
When we assume random lensing alignments, the distribu- 
tion of the CSL values depend only on the properties of the 
lensing profiles. For a single universal cusp slope value, the 
numerical formulation for the CSL probability is presented 
in section [2~2l 



2.1 Condition for triple image lensing 

An axially symmetric lens potential produces visible triple 
image lens when images at ki and k$ are sufficiently close 
to the radial critical curve at fc cr i = v. Otherwise, the image 
kz inside the caustic curve becomes unobservable because of 
strong demagnification, while image k\ remains magnified, 
and a double image lens system is observed. 

Images nearby the critical radius k crl are produced by 
source images that are located near I ~ Z ma x, where Z max 
is the maximum source coordinate for producing multiple 
images (strong lensing). Thus, the triple image condition 
demarcates out part of the strong lensing region. In other 
words, for double image lensing, the source coordinate I must 

have I < /triple < ^max- 

The condition for triple image lensing can be derived 
by examining the magnification ratio /ii//i3 of the images 
ki and £3. This is done by expanding the lens equation at 
k — v, setting the condition /ii//i3 < r and solving it for 
the corresponding source coordinate. Here r is the limiting 
value for the observable magnification ratio. 
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From the lensing theory presented in 
iMutka fc Mahonen I ||2006h . the maximum source coor- 
dinate has value 



,1 



a + 2 



(16) 



where v = (a + 2)~ 1// ' a+1 - ) and a is the cusp slope of the 
density profile in equation ([1]). 

We normalize the lens equation 

l -=k(l-\k\ a+1 ) (17) 

with expression (I16p in case of— 2^a<— 1. Now the lens 
equation can be expressed in new coordinates y = Z/(gL max ) 
and x — k/v as 



y = 



a + l 



(a + 2-x a+1 ). 



(18) 



We continue our derivation by expanding the re-normalized 
lens equation at k = k crl = v 4$ x = 1, which gives us 



V 



1 a + 2 > i^s 
1 7T- (a:- 1 ) • 



(19) 



By solving this approximation we get the image location 



x = 1± 



a + 2 



(20) 



from a source nearby the critical curve. Here the positive 
sign corresponds to image ki an d negative to k^. Aga i n, from 
the lensing theory presented in IMutka fc Mahonen I {2006), 
the magnification ratio condition for the triple image lensing 
can be expressed as 



Mi 

M3 



\y-xi\ 



< T. 



(21) 



By inserting solutions for xi and X3 from equation (1201) to 
this expression and solving for the source coordinate y, the 
condition for observable triple image lensing becomes 



y > 1 



a+2 (T- 1 
r + 1 



(22) 



The power of two factor in parenthesis that depends on r, 
approaches very quickly to the unit value if demanding ten 
or a hundredfold magnification ratio for the limit of observ- 
ability. Because determining the exact value for the mag- 
nification ratio limit is impossible, we set this factor equal 
to one. Thus, practically without any loss of accuracy, the 
source image condition for the triple image lensing can be 
generalized to 



a I av a+l triple 
V > ~ 2 ^ q > ~ a + 2 ~ q 



(23) 



As before, we assume that the source images are uni- 
formly distributed at the source plane inside radius Z/ max . 
Now the source image coordinate limit (|23[) leads directly to 
a simple graphical probability for triple image lensing 



Ptriple(a) 



a 

T 



(24) 



that depends only on the cusp slope parameter of the den- 
sity profile {TJ. The expression (|24[1 for probability of the 
triple image lensing can be considered also valid with reason- 
able values of ellipticities for triaxial haloes. Corresponding 



treatment can be made for haloes with elliptic isocontours 
following the profile ((T}. Note, that the probability (|24[) ap- 
plies only on proportional share of triple image lenses among 
already identified lens systems, not on arbitrary galaxy in 
the Universe. This is why the q term dependency vanishes 
from the equation. Additionally, this result does not account 
for the finite angular resolution in the real world observa- 
tions. 

The remarkably simple condition (|23p and probabil- 
ity (|24[) have direct consequences on the general proper- 
ties of the lensing potential and properties of dark matter 
haloes. According to the expression (|24p . double image lens- 
ing and triple image lensing have equal probability when 
a = a cq = — y2 ~ —1.4. If the cusp slope value a > a oq , 
triple image lenses are more common than double image 
lenses. The fact that there are only three known cases with 
visible inner caustic image out of roughly hundred impli- 
cates that the cusp slope in the dark matter profile must be 
closer to isothermal value a ~ — 2 than a cq . 

In fact, if we set P tr i P io = N tIip i<./N — 0.03 and solve 
the cusp slope parameter from the equation (|24[) . we get 
a w —1.97 - a result rather close to the quoted value for the 
population HI in this paper. 

2.2 Theoretical CSL probability 

Here we show how to calculate a numerical probability func- 
tion for the CSL limit values. We assume that the source im- 
ages are uniformly distributed at the double image region of 
the source plane limited by the equation (|23j) . i.e. the radial 
source coordinate 

"=5£: < -? <1 - (25) 

We assume a single population of lenses with a universal 
cusp slope value a. All the slope parameters used in the 
following calculations have —2 < a < — 1. Special cases 
a — — 1 and a = —2 could be treated in similar way, but it 
is not necessary at this point. 

The differential probability for having a lens with value 
a C sL is corresponding to an area of a ring with an inner 
radius y and an outer radius y + dy: 



dP : 



8^dy, 



(26) 



where y = j/(acsL;a). This equation can be integrated in 
order to estimate the probability of a double image lens to 
have a2 < a C sL < oti: 



P(a 2 < a C sL < ai) = 4 



y 2 (a 2 ;a) - y 2 (a 1 ;a) 



(27) 



Although it could be possible to write an approximation 
for the function y = y(a C sL;a) or rewrite the differential 
(1261) using other variables, we calculate it here numerically. 
For this purpose we employ the magnification ratio M = 
M(a, 9) of the images k\ and k% produced by the source 
image at as presented in equation (|13[) . Here 9 = fci/fo is 
the image coordinate ratio. 

We start calculating the y(acsL;a) for each Ocsl by 
solving the equation 

$1 
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numerically for #csl- The acquired value is subsequently in- 
serted to 



(29) 



M(a C sL,0csL) = M(a, 



which is then solved numerically for 0. Equation (|29[1 has 
two solutions 6 — 0( u )(o!csl) and 8 = 9^) (ckcsl)- When 
these solutions are inserted into the lens equation 



V = y{9,a) 



i + e a + 2 



l/(t*+l) a + l 



+ 1' 



(30) 



source coordinate values ?/(„) (o;csl; a) and j/(d)(acsL; a) are 
acquired. Because we must limit the solutions to the double 
image cases, resulting yr u -\ or values are not allowed to 
exceed the limit (12311 , 

With these two solutions, the probability becomes 



P(a 2 < cucsl < ai) 
y(d)(«i; a ) + y(«)("2;a; 



(31) 



from radio and optical observations were used to calculate 
mean flux ratio weighed with measurement errors in obser- 
vations (columns four and five). The flux ratio is calculated 
by dividing the closer (dimmer) image flux with the further 
(brighter) image flux. Columns six and seven contain eval- 
uation for redshifts of the lens candidate (if present) and 
the source correspondingly. The eighth and ninth columns 
in the table contain radio and optical cusp slope limit values 
(Qcsl). Note, that these may have different values because 
the radio and optical sources have different locations at the 
source plane. The last column lists the relevant references 
pointing to the source of the data and/or the publication 
reporting the discovery of the lens system. 

After the preliminary Monte Carlo testing, we will use 
the calculated CSL values in the final statistical analysis. 
Because the optical data is more abundant, we have chosen 
to prefer optical values, and use radio measurements only if 
the optical data is not available. 



3 DOUBLE IMAGE LENS SYSTEMS 

For this study, we composed a catalogue of double image 
lens systems, which have a configuration resembling axially 
symmetric lensing. The catalogue is based on JVAS-CLASS 
radio survey and optical CASTLES survey and it contains 
44 lens systems, see table [1] and relevant references. Four of 
these systems had no data available during the conducted 
background research, which reduces number of useful sys- 
tems down to 40. 

The Jodrell Ba n k - VLA Astro-metr i c Survey (JVAS, 
Patnaik et all 1 19921 ; iBrowne et all Il998l : IWilkinson et all 
19981) and the Cosmic Lens Al l Sky Survey (CLASS, 



Jacks on et al ] 1 19951 : iMvers et al.l Il995l ) were originally 
geared towards finding new gravitationally lensed systems 
in order to determine the value of the Hubble's constant 
Ho- During the surveys, they observed roughly 10000 flat 
spectrum radio sources. At the moment these surveys has 
detected 20 strong lens systems. The Cfa-Arizona Space 
Telescope LEns Survey is an ongoing HST survey of all the 
known strong lenses and lens candidates. The main goal of 
the survey is to create a n uniform high qua lity photometric 
sample of these objects (jMufioz et al.lll998l m 

Our subsample of these surveys was chosen according to 
the following criteria: the lens system must have two clearly 
separate lensed images corresponding to a quasar source. 
Triple and quad lenses, Einstein rings and clearly non- 
axially-symmetric cases were excluded. Additional lensed 
components, such as radio jets were allowed as long as the 
lensed point source was present. 

Table [T] summarizes properties of the accepted lens sys- 
tems. After the name of the lens system, absence of a lens 
object candidate is indicated by a '-' in the second column of 
the table. If there are additional components at the imme- 
diate surroundings of the lens system, it is indicated with a 
'x' in the third column. All the available flux measurements 



1 C.S. Kochanek, E.E. Falco, C. Impey, J. Lehar, B. 
McLeod, H.-W. Rix, the CASTLES survey has a website 
http : //www . cf a. harvard. edu/castles/ 



4 MONTE CARLO TESTING 

In order to test the CSL-analysis method against different 
perturbations and investigate the error sources, we have 
developed a computer code that generates mock lens cat- 
alogues. Our code uses rejection method with base distri- 
butions that are randomly sampled for initial redshifts, lens 
masses and source luminosities. Sampled values are fed to 
the lens equation solver that accounts for the lens asymme- 
tries and other more sophisticated properties of the lensing 
event. Generated lens system is accepted or rejected accord- 
ing to the observed configuration, selected magnitude limit 
and minimum image separation threshold. 

We use differential optical depth 4j (Ml, z s ), as a func- 
tion of lens mass M L and source redshift Zs for the base 
lens distribution (|Ofek. Rix fc Maozl 120031 ). Our generator 
uses similar morphological mix for the lens galaxies, and 
we have chosen parameter values U = —0.20, P = 1.20 and 
/ E = 0.95 for the distribution. See lOfek. Rix fc Maozl (|2003h 
for further details, and information on other parameters af- 
fecting the distribution. 

The source objects are sampled from a double power 
la w luminosity function < fr(L, z s ) of quasars as presented 
bv lWvithe fc Loebl (|2002l) . They use a luminosity function 
with pure luminosity evo lution and additional br eak for very 
high redshift objects. See lWvithe fc Loeb I (|2002l ) for further 
details. 

Our parameters generate reasonably good samples of 
lens systems, as can be seen from figure [4] We compare our 
mock data against observed redshift distribution of lenses, 
redshift distribution of sources, image separations and lens 
masses inferred from the Faber-Jackson relation. We use 
Kolmogorov-Smirnov two-sample test providing a probabil- 
ity that the samples are drawn from the same background 
distribution. When observed lens object magnitudes are con- 
verted to absolute luminosities for mass estimation, evolu- 
tionary a nd K-corre c tions are a pplied to the data as pre- 
sented by Poggiantij (|l997l ) and iBicker et al. | l|2004l ). Note, 
that the size of the observational sample is changing in each 
panel. The reason for this is that we use all the double image 
lens data available for testing, and some lens systems have 
missing data. 
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Figure 4. Lens data generated by the simplest model in our Monte Carlo code. Panel a) presents distribution of the lens redshifts. The 
sampled redshift distribution contains some lenses with redshifts z L > 1 (dashed line), but they vanish if the lens object magnitudes are 
accounted for (solid line). Shaded region outlines 2cr-region corresponding the Poisson noise proportional to the observed lens sample size. 
Over-plotted histogram describes redshift distribution of observed lenses. Panel b) shows similar redshift distribution for the sources. 
The image separations in arcseconds (9) are presented in panel c) and sampled lens massed compared to the lens masses inferred from 
observed magnitudes for the lens objects are plotted in panel d). Note, that here the symbol 9 for the image separation should not be 
confused with the image coordinate ratio in equation iY.it . The observed data is from the CASTLES survey catalogue. Each panel lists 
a Kolmogorov - Smirnov test probability P against observed data, ratio of the bins with more than five data points (ris) to the total 
number of data points (n to t), binsizes (Azl, Azs, A#, AMl), and total number of lens systems (n to t). 



In figure[4] the panel c) presenting image separation dis- 
tribution has low Kolmogorov - Smirnov probability. Note, 
that here 9 is image separation in arcseconds, that should 
not be confused with the image coordinate ratio in equa- 
tion 1131 The mock data distribution follows the observed 
image separations when 9 > 1", while lenses with tighter 
image separations are over represented. We believe that this 
discrepancy is a result from an observational bias. Because 
the observed lens systems are not from a single systematic 
survey, the lenses with small image separations are under- 



represented in observations (lenses with image separation 
> 1" are easier to detect). It is clear, that the lens data ac- 
quired with a large array of different instruments is far from 
being statistically uniform. Therefore, chosen values for the 
magnitude limit and image separation limit were determined 
from fitting to the data. 

Although we conducted a parameter sweep for mag- 
nitude limit Mum) image separation limit 9 Um and param- 
eters U, P, / E , it should be emphasized that we do not 
try to find the best fit for the base-distribution parameters. 



10 P. T. Mutka 



Table 2. Numerical values for the parameters used by our Monte 
Carlo code. The upper part of the table lists parameters for the 
Faber-Jackson relation and redshift evolution in the differential 
optical depth for lensing d"r/dz(M L , z s ). The lower part contains 
parameters for the quasar luminosity function <&(L,2 B ). See text 
for the quoted references on the further details about the param- 
eters. 
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We use these values only for providing as realistic test data 
for the CSL-analysis method as possible. It should also be 
noted that even though the generated data is not exactly 
fitting the observed sample, it should contain similar inter- 
nal correlations as the real data. Therefore it can provide 
decent test cases for our method. We use bolometric mag- 
nitude limit is Mii m = 30.5 mag and image separation limit 
#um = 0.2". Table [5] summarizes chosen parameters for the 
base distributions. For further e xplanation about physica l 
mea ning of the parameters , see lOfek. Rix fc Maoa (|2003l ) 
and lWvithe fc Loeb I l|2002j ). 

In the following subsections, we explore the effects from 
the ellipticity of the lensing potential, the substructure 
within the lensing halo, and the source variability coupled 
to the time-delay that can distort the magnification ratios 
from the lensing events. This is done by creating mock lens 
catalogues with different degrees of perturbations and per- 
forming the CSL-analysis on them. 

4.1 The Basic Concept 

The concept of the CSL value is to disentangle cosmology 
and geometrical configuration of the each lens event from the 
properties of the lens halo profile. These factors cancel out in 
the derivation of the CSL equations (|14[) and (|15|l . However, 
they do have an indirect influence on the underlying sample 
of the lens systems and can affect the statistics of the CSL 
values. 

The disentanglement of the lens halo properties from 
the other factors in the CSL analysis is possible because 
we assume simple lensing properties - the lens potential is 
axially symmetric and the source images are uniformly dis- 
tributed on the strong lensing region of the source plane. If 
there is a strong bias in the actual source image locations on 
the source plane, it can distort the resulting CSL distribu- 



tion. Such a bias can emerge in a magnitude limited sample, 
where sources at the high magnification region of the source 
plane are over-represented in the total sample. 

We have constructed a numerical probability function 
for the CSL values (oicsl) with an universal cusp slope a 
(see appendix 12. 2p . We fit this function to a histogram of 
Qcsl values calculated from the (mock) data. In the fitting 
procedure, we use the Levenberg-Marquardt method and 
standard ^-statistics. The fitted parameters are the uni- 
versal cusp slope value a and normalization parameter n. 
Figure [5] illustrates such a fit. 

The universal density profile has been one of the main 
assumptions in studying the properties of dark matter 
haloes. A population of lenses with a single cusp slope pa- 
rameter a produces a distinctive CSL-value distribution (fig- 
ure [5J that is very different from completely random pro- 
files. Figure [6] shows CSL distribution of a lens population 
with uniformly distributed random profiles with a in range 
[—2,-1] (hatched area), and second lens sample with pro- 
file slope parameters in same range with normal distribution 
(expectation value (a) = —1.5 and a a — 0.15, solid line). 

When performing fits on bin ned data, there are some 
issues to be considered, see e.g. iHumphrev. Liu fc Buotel 
(2009). We tested several goodness-of-fit functions and bin- 
sizes in order to find such a configuration that can re- 
cover the universal cusp slope value from the data as of- 
ten as possible. We have chosen to use the previously men- 
tioned goodness-of-fit function and bin the data in 14 bins 
in range o C sl = [—2.0, —1.07]. By using our binning selec- 
tion and goodness-of-fit function, we were able to recover 
the universal cusp slope a with 2a accuracy in most of the 
cases (cusp slope value a = [—2.0, —1.1] and catalogue size 
Ni ms = [30,100]). 

The quoted error limits for the fitted parameters and 
corresponding rejection probabilities are calculated with the 
Monte-Carlo method. We use a bootstrap method to create 
hundred data sets that are fitted using the same procedure 
as the original data. These fits are employed for error es- 
timates and rejection probability calculation. This is a well 
known procedure that has been e mplo yed b y several authors, 
see e.g . papers by iBuote et al. I (|2003l ) and lHumphrev et al.l 
(2006) for further information. 

There are certainly more suitable values for binning that 
can be adapted for certain sizes of samples and values for 
the universal cusp slope a. However, we emphasize a good 
overall performance with wide range of sizes for the lens 
samples and parameters. If not mentioned otherwise, we use 
these choices in the analysis presented in the following sub- 
sections. 

4.2 Lens Ellipticities 

The CSL-analysis is based on an axially symmetric lens 
model. In reality, the observations and the N-body simu- 
lations indicate t hat the dark matter haloes are triaxial el- 
lipsoids, see e.g. lHavashi. Navarro fc Springel I (|2007f ) and 
references therein. Here we examine the effects from asym- 
metries that can distort the CSL-analysis results by chang- 
ing the observed magnification ratios. 

We have generalized our axially symmetric lens model 
to the elliptic lens potential case. The isopotential surfaces 
of the elliptic lens follow similar surface density as the ax- 
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Figure 5. Fit of a re-normalized theoretical probability function 
for the CSL distribution of a lens population with cusp slope 
a = —1.8. Jagged histogram and corresponding diamond symbols 
with error bars present histogram data for ocsl values calculated 
from the mock sample. Error bars correspond to 0.1 magnitude 
absolute error in photometry. Thick solid curve represents fitted 
CSL-distribution, with a (dashed line) and 2a (shaded region) 
outlining the Poisson shot noise error region. The fitted cusp slope 
value for the data is a = —1.82 ± 0.01 and n = 6.9 ± 0.3 with 
P < 0.01 rejection probability. 



-2.0 -1.8 -1.6 -1.4 -1.2 -1.0 
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Figure 6. If the cusp slopes of lensing haloes are random, char- 
acteristic CSL-distribution (see figure [5} disappears. This figure 
illustrates how a population of lenses with random cusp slope 
values appears in the CSL histogram. The hatched area is a his- 
togram of a lens population with cusp slope values that have uni- 
form distribution. Thick line presents histogram of a lenses with 
normal cusp slope values with an expectation value at (a) = —1.5, 
and variance a a = 0.15. Error bars represent 0.1 magnitude ab- 
solute error in photometry. 



ially symmetric lens model, equation Jl]). Our elliptic lens 
model approaches the axial model with eccentricity value 
e — > 0. The lens equations were deriv ed from numerica l in- 
tegrals for elliptic lensing present ed in ISchramml (|l990l ) and 
ISchneider. Ehlers fc Falcof l| 19921 ). 

The elliptic lenses can produce four image lenses, when 
the source image is located inside the inner diamond shaped 
caustic curve. This configuration is related to the Einstein 
ring in the case of axial symmetry. Our Monte Carlo code 
accepts only double image lens systems into the sample, i.e. 
all the source images must be outside the inner caustic. The 
size of the area at the source plane demarcated by the inner 
caustic curve depends on the cusp slope a of the profile 
and on the eccentricity e of the halo. In general, a larger 
inner caustic area means stronger deviations from the axially 
symmetric lensing. 

Strongly cusped lens profiles can tolerate high degrees 
of eccentricities without significant difference in the image 
flux ratio to the axially symmetric case. Deviations from 
the axially symmetric model grow progressively larger with 
shallower profiles and increasing eccentricities. This trend is 
illustrated by figure[7] in which we have examined an unreal- 
istic scenario where all the lenses have constant eccentricity. 
We have constructed lens catalogues with cusp slope values 
a = [—1.9, —1.75, —1.6, —1.45, —1.3] and varied eccentric- 
ity of the lensing potential. We investigated eccentricities in 
range e = [0, 0.6] with 0.05 increments. 

In figure [71 each cusp slope value is represented with 



a horizontal dashed line. Plotted data points correspond to 
fits to mock catalogues including error estimates with each 
value of eccentricity. Data points are plotted until the el- 
lipticity makes the CSL distribution unrecognisable (cutoff 
vanishes and tail becomes erratic). The figure clearly shows 
the aforementioned trend for cusped profiles; the fit tends 
to increasingly overestimate the cuspiness of the profile with 
increasing halo ellipticity. The cutoff comes earlier with shal- 
lower profiles. 

We model the lens potential eccentricities with log- 
normal random distribution with expectation value (e) and 
va riance <j e - As a reference case, we follow the treatment 
by Irluterer. Keeton fc Ma I ([2005) and adopt the distribu- 
tion of ellipticities measured fo r 379 early-type galaxies by 
|j0rgensen. Franx fc Kiaergaardl l|l995h . The distribution has 
mean (e) = 0.31 and a e — 0.18, with upper limit e < 0.6. 

Our results indicate that, regardless of realistically el- 
liptical lens haloes, the modelled CSL value distribution can 
recover the universal cusp slope up to slope value a ~ —1.7 
and give a rough guess for the exact value up to a ~ —1.5 
that is always an overestimate. See figure [5] 

4.3 Substructure 

The dark matter haloes should have several sub-haloes that 
generate irregularities in the lensing potential. These irreg- 
ularities can corrupt the macro-lens magnification ratio and 
distort the CSL-distri bution. Here we ha ve chosen to follow 
similar approach as in lRozo et al.l l|2006l ). We model the ef- 
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Figure 7. We have compiled lens catalogues with cusp slope val- 
ues a = [—1.9, —1.75, —1.6, —1.45, —1.3] with constant eccentric- 
ities. The lens catalogues have 100 lens systems. The x-axis mea- 
sures eccentricity e and the y-axis cusp slope a. Each data point 
corresponds to a fitted CSL distribution with error estimates. All 
the fits with recognisable cutoff in the CSL value distribution 
are plotted. If the eccentricity is further increased, characteristic 
shape (figure [5J of the distribution vanishes. 



fects from the substructure on the magnification of images 
with iV linearized perturbers contributing with Sfu to the 
total perturbation 5. When contribution from all the per- 
turbers is added together, we get the total perturbation 



5 = 



2\fi\ 



(1 — k) 2_J 5Ki — 7 YJ S*fi cos(2<j>i) 



(32) 



where k and 7 are macro-lens convergence and shear, 
whereas is macro-lens magnification. Each perturber is 
located at radial coordinates 9i and (pi, and convergence per- 
turbations and shear perturbations 5fi are correspond- 
ing contributions from the ith perturber. It is assumed that 
\5ni\ <C 1, thus the astrometric perturbations by the sub- 
structure are negligible and only the changes in the image 

fluxes are considered. 

We follow the analysis bv lDalal fc Koch anck ( 2002]) and 
iRozo et al.l (|2006l ) by employing a pseudo-Jaffe density pro- 
file 



p(r) oar (r + a 2 ) 



(33) 



for individual perturbers. Here r is radial coordinate and a i s 
an effective tidal radius. As presented in lRozo et al. I i|2006l ). 
we parametrize the tidal radius with X a as 



(34) 



where b is the Einstein radius of the perturber and 6 H the 
corresponding radius for the macro-lens model. Convergence 
and shear contribution (5/tj and Syi) from the ith perturber 
at polar coordinates (9i — kiX B r B / D^, <j>i) is 



Figure 8. Cusp slopes from the fitted CSL values using log- 
normal random distribution for the lens halo eccentricities (see 
text). The x-axis has the cusp slope value used in the mock cat- 
alogue generator, and the y-axis has fitted cusp slope value. Di- 
agonal dashed line marks the exact fit. Each data point denotes 
a result from a CSL distribution fit to a compiled lens catalogue, 
with error estimates. Lower panes shows fit residuals. Fit gets pro- 
gressively worse with more shallow profile until at a ~ —1.45 the 
characteristic shape of the CSL value distribution is too distorted 
in order to make a sensible fit to the data. 



Sni — b 
and 

5ji = b 



1 1 
20~ _ 2$ 



1 



1 



a a 



29i + 2£ 9° ' ' 



(35) 



(36) 



The mass of the perturber is m = -KabT, CI D^, while b is 
defined by 



b „ a 
b= 1+ b- 
and £ is 







1 + 







(37) 



+ a 2 . (38) 

When a — > 00 the pseudo-Jaffe lens profile approaches SIS 
profile, and when a — > the profile corresponds to a point 
mass lens. 

Mass spectrum s of the perturbers is assume d to follow 
ds/dm oc m 13 , where we choose /3 = —1.8 as in iGao et alj 
(2004). The amount of substructure is defined by a ratio 
of the substructure surface density to the critical density 
/ sub = 2E a /E cr . The mass spectrum is parametrized also by 
the cutoff at the substructure mass value m max that is chosen 
to be 1 per cent of the mass within the Einstein radius (mass 
M E ) of the lens object. 

We account the substructure perturbations to the im- 
age magnification in the synthetic lens catalogue generation 
as follows: the magnification change of the each image is es- 
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Figure 9. Cusp slopes from the fitted CSL values including sub- 
structure perturbations to then magnification. The x-axis has the 
cusp slope value used in the mock catalogue generator, and y-axis 
has fitted cups slope value. Diagonal dashed line marks the exact 
fit. Each data point denotes a result from a CSL distribution fit 
to a compiled lens catalogue, with error estimates. Lower panes 
shows fit residuals. 



timated by randomly generating 1000 perturbers uniformly 
distributed at the lens surface with a surface density that is 
determined by the / aub . The total perturbation Stot is then 
estimated from the equation (|32p with negative surface mass 
component in order to guarantee the mass conservation 



5 t ot = S - - «;)/„ 



(39) 



The perturber masses follow the previously discussed power 
law mass spectrum with index j3. The total perturbation 
Stot is then added to the total magnification of the image, 
corresponding to the macro l ens shear an d conv ergence. 

Again, we follow the Roz o et all ([2006) with the 
parametrization of the substructure and choose tidal ra- 
dius parameter X a = 4.0, and substructure surface density 
/ sub = 0.005. 

We study the goodness of the CSL fit as a function 
of the universal cusp slope of the generated lens catalogue 
with constant strength for the substructure. The lens cata- 
logues are generated with an axially symmetric lens model 
with cusp slopes a ranging in [—1.95,-1.45]. The results 
are summarized in figure [9] All the fits are very close to the 
original value, although it would seem that fits get slightly 
worse at higher values for the cusp slope (a > —1.7). 

Although the substructure can significantly alter the 
observed flux ratio of the images, this can happen only at 
the vicinity of the critical curves for lensing (at the source 
coordinate I ~ or I ~ £ m a, x )- These kinds of lens systems 
produce either strongly magnified Einstein rings or triple 
image lenses that are excluded from the analysis by the se- 
lection criteria. Furthermore, because such lens systems have 
incidence probability proportional to the square of the ratio 



1/ltnxK, they should be comparatively rare when the source 
images are uniformly distributed at the source plane. 

Here it should be noted, that a similar treatment could 
be applied for micro-lensing. However, because the CSL- 
limit analysis is based on double image lenses where the im- 
ages are generally not close to the critical curves, we assume 
that the effects from micro-lensing can be safely neglected. 

4.4 Time-delay 

In the double image lenses, the lensed images represent the 
source at different moments in time in the rest frame of the 
source. This is a consequence of a different light propagation 
distance for each image (geometric time-delay) and different 
gravitational potential along the light ray trajectory (poten- 
tial time delay). As a consequence, the observed flux ratio 
is different from the magnification ratio of the images if the 
source flux is varying at suitable timescale. This can distort 
the acquired CSL-distribution. 

We calculate the time delay from the lensing po- 
tential generated by our parametrized density pro- 
file. For definitions of time delay, see for example 
ISchneider. Ehlers fc Falcol l| 19921 ). The source variability is 
modelled by generation power law noise, that is parameter- 
ized according to the obs ervational relations published in 
IVanden Berk et al. I (|2004h - 

Source flux variations are modelled by generating power 
law variations Am around the baseline magnitude of the 
source sampled from the quasar luminosity function. The 
power spectrum of the magnitude variations is 



(40) 



where ui is frequency, r\ power law index and $o correspond- 
in g normalization. On gener ating power law noise, see paper 
bv lTimmer fc Konid i| 19951 ). 

Noise with characteristics defined by (|40p have variance 



^Am = 2 J dLU<f>(iu) 

2vr/ 



2$ (27r)-" !_„ 



(41) 



where / is the minimum frequency to be considered. Corre- 
sponding covariance function for the magnitude m = m(r) 
is 



B(At) = (m(T)m(T + At)} = {Am z ) 
= I -!o;cos^r)$(^) = i$or''- 1 A(» 7 ,/), (42) 



OO 

/ 

2ir/ 

that is written at time r with time lag At. In previous 
equation we define A(r], f) as 

A(f],f) = 

e- i ' r(1 -' ?)/2 r(l - n, i2nf) + e -( 1 -")/ 2 r(l - v , -i2wf), (43) 

where T{x,y) is the incomplete gamma-function. The vari- 
ability function V = V(r) and model for it are defined as 



V(T) = x l^(Am)*-(** /n )=V 



(44) 



where Vo and to, that are determined from the observations, 
parametrize the model. Here (of/a) i s the signal to noise ratio 
of the variations. 
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When (Am 2 ) ~ 7r(Am} 2 /2, we get from the definition 
of the variance and equations (I41|l and (|42|) the relation be- 
tween power spectrum parameters and the variability func- 
tion. The power law indices are related as 

7? = 1 + 2 7 (45) 

and normalization as 

^(»?>/)V 

Here it should also hold (cr 2 /n ) = 7r | o^m I • 

The normalization $o and the power law index rj charac- 
terizing the power law no ise are calculated from the v ariabil- 
ity function presented bv lVanden Berk et al. I (|2004h . where 
they study photometric variability of roughly 25 000 quasars 
in the SDSS survey. 

In their paper, they estimate the variability function 
power law index as 7 = 0.246 ± 0.008 and characteristic 
time scale to = (5.36 ± 1.46) x 10 5 d. The normalization of 
the variability function (|44[) is parametrized as 



source variability (N, ens =100) 



V = v(M)v(X K )v{z) 



(47) 



with the absolute magnitude of the quasar M, observed rest- 
frame wavelength A R and redshift of the quasar z. Here we 
set 



v[M) = 10 W2 - S 

in which /3 = 0.246 ± 0.005. Correspondingly 

v(\ R ) = a e~ x/Xa + ai 



(48) 



(49) 



with a = 0.616 ± 0.056, ai = 0.164 ± 0.003 and Ao = 
988 ± 60A. And finally 



v(z) = (0.019 ± 0.002)2 + (0.037 ± 0.005). 



(50) 



The time lag At is acquired from the time delay equa- 
tions, and we choose rest-frame wavelength corresponding to 
the observed B-band wavelength. The absolute magnitude 
is acquired from the sampled luminosity of the quasar. 

We convert the parametrized variability function to the 
power law noise index r\ and corresponding normalization 
$0, and obtain random change in the flux of the source by 
generating similar power law noise. This changed flux (that 
equals the baseline flux sampled from the quasar luminosity 
function with the random change for the each image) is then 
correspondingly magnified by the lens and inserted into the 
standard mock lens catalogue generator. 

We study the effects from the time delay coupled to 
the source variability by generating several mock cata- 
logues with an universal cusp slope values a ranging in 
[—1.95, —1.45]. The mock catalogues are generated with an 
axially symmetric model. The fitted CSL values are pre- 
sented in figure [TD] 

Inclusion of the source variability seems to slightly 
strengthen the acquired cusp from the lensing data. Al- 
though all the fits are reasonably good, strongly cusped cat- 
alogues are most affected by the source variability. 



4.5 Combined error sources 

In the previous subsections, we have examined the strength 
of several distortions on the CSL statistics with the simplest 
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Figure 10. Cusp slopes from the fitted CSL values including 
source variability induced time delay perturbations to the image 
flux. The x-axis has the cusp slope value used in the mock cata- 
logue generator, and y-axis has fitted cusp slope value. Diagonal 
dashed line marks the exact fit. Each data point denotes a result 
from a CSL distribution fit to a compiled lens catalogue, including 
error estimates. Lower panel shows fit residuals. 



possible model. Here we add all the contributing factors (el- 
lipticities, halo substructure and time-delay) together in or- 
der to see what is the applicability of our method. 

We study the goodness of the CSL fit as a function 
of the universal cusp slope value in the mock catalogues. 
We use elliptical lens model with random log-normal ec- 
centricities as in subsection 14.21 The modelled lens systems 
possess similarly parametrized substructure as in 14.31 and 
source variability corresponding to the previous subsection 
4741 

The results of the fits are presented in figure [TT1 It would 
seem that a simple axially symmetric model with an assump- 
tion of uniform source image distribution at the source plane 
can capture essential properties of the double image lenses 
in the presence of higher order effects when calculating the 
CSL statistics. 

In general,the CSL goodness-of-fit x 2 -surface contains a 
local minimum inside each bin. In the fits that are presented 
here, we have always chosen the global minimum. However, 
in some cases there are two local minima that have almost 
the same % 2 values resulting two equally good fits in sta- 
tistical sense. Such problems can be resolved by manually 
adjusting the binning to cover the range of cvcsl values in 
an optimal way. We have chosen to use a fixed binning that 
produces fits that should become progressively poorer when 
the universal cusp slope approaches a = — 1 because the 
data is assigned to fewer bins correspondingly. However, the 
leftover zero bins are fitted into the model with estimated 
error bars, thus the recovered error bars for a do not change 
significantly. This is a conscious choice that was made in 
order to study the error sources as a function of the cusp 
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Figure 11. Cusp slopes from the fitted CSL values with all the 
examined perturbation sources (lens ellipticity, substructure and 
time delay effects coupled to the source variability). The x-axis 
has the cusp slope value used in the mock catalogue generator, 
and y-axis has fitted cusp slope value. Diagonal dashed line marks 
the exact fit. Each data point denotes a result from a CSL distri- 
bution fit to a compiled lens catalogue, including error estimates. 
Lower panel shows fit residuals. 

slope in a consistent way. In principle, an adaptive binning 
should produce better results, but as the real data is spread 
over rather wide range of the cvcsl values, we hold to our 
previously adopted binning scheme. 



5.1 Single population model 

We start our analysis with an assumption of a single univer- 
sal halo profile. The acquired cusp slope value for this lens 
population is a = —1.95 ±0.02. The value of the normaliza- 
tion parameter is not interesting, because it is simply scaling 
the probability function to the size of the sample. The re- 
sulting fit is presented in figure [T2l The acquired value for 
the cusp is close to the value of an isothermal sphere (SIS). 

All the data points below «csl ~ — 1.5 fall well into the 
one Poisson-cr range of the model. However, at cvcsl > — 1.5 
there are two data points well outside of the 2cr-Poisson noise 
range. These data points represent roughly 18 per cent of the 
total sample of lens systems. All corresponding lens systems 
cannot belong to the same profile group as the acquired 
fit because perturbations have a tendency to steepen the 
acquired cusp slope. One explanation could be that these are 
strongly sheared systems, but there are several arguments 
against this possibility. 

The fact that the deviating lens systems are rather 
tightly clustered around value a ~ —1.5 does not support 
sheared lens scenario. Randomly sheared systems should 
produce wider distribution, with similar characteristics as 
random cusp slope values (see figure [5} . Strongly cusped 
haloes are also more resistant against the deviations from 
the axial symmetry than the haloes with shallower cusps 
(see section I4.2|l . Additionally, sheared systems that have 
strongly deviating magnification ratio from the axially sym- 
metric value, very often produce more than two images in 
similar fashion as strongly elliptic haloes. Such lens systems 
are excluded from the analysis by the selection criteria used 
to construct our lens sample. 

A single universal cusp slope value cannot fully explain 
the observed CSL-value distribution. Therefore, we investi- 
gate also a possibility that this deviation is a signature of 
a second population of haloes with their own characteristic 
density profile. 



5 FITTED CSL STATISTICS 

In this section, we apply the statistical analysis of the CSL 
values to the real double image lens data. We use our dou- 
ble image lens sample described in section [3] to study the 
distribution of the CSL values. The calculated CSL values 
for each lens system are presented in the table [T] These val- 
ues are fitted to the theoretical distribution (see appendix 
I2.2|l . The fit produces estimates for an universal cusp slope 
value a and corresponding normalization n proportional to 
the size of the sample. 

In the fitting procedure, the error limits for the his- 
togram were created using the Monte-Carlo method. Thou- 
sand distributions were constructed with simulated Gaus- 
sian errors according to the measurement errors in the flux 
ratios. The data points in the presented histogram are mean 
values, and the corresponding error limits are the standard 
deviation calculated from these constructed histograms. In 
some cases values at the end points of the distribution had 
zero error limits. Those were estimated upwards with mean 
of non zero errors within the distribution. We use the same 
setup for the data as for the Monte Carlo testing described 
in the previous section. The details are presented in section 

ED 



5.2 Dual population model 

The dual population model (populations are abbreviated as 
HI and H2) is implemented by adding two individually nor- 
malized distributions together, producing a model with four 
parameters. This model is fitted to the observational data 
in a similar fashion as the single population model in the 
previous subsection. The fit produces estimates for the uni- 
versal cusp slope values a H i and a H 2 for both populations, 
and the corresponding normalization values n H i and n H 2- 

The detailed fit to this data is presented in figure 1131 
The acquired values are cusp slope a H i = —1.95 ± 0.02 and 
normalization n H i = 2.5 ± 0.7 for the population HI and 
for the population H2 correspondingly a H2 = —1.49 ± 0.09, 
nn2 = 0.5 ±0.3. Normalization values n m and n H 2 gives the 
relative abundance of population HI and population H2 lens 
objects in the sample. Our results indicate that 83 per cent 
of the lens objects belong to the population HI and 17 per 
cent to the population H2 with 4 per cent error margins in 
both values. This suggests that on average roughly one sixth 
of the lens objects capable of strong lensing belong to the 
second population with shallower cusp slope. 

The dual population fit is better than the single pop- 
ulation fit, which is to be expected because number of free 
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Figure 12. A single population fit of a re-normalized theoreti- 
cal probability function for the CSL value distribution from an 
observational lens catalogue. Jagged histogram and correspond- 
ing diamond symbols with error bars present histogram data for 
ocsl values calculated from the data. Error bars correspond to 
measurement errors in photometry. Thick solid curve represents 
fitted CSL-distribution, with <r (dashed line) and 2a (shaded re- 
gion) outlining the Poisson shot noise error region. The fitted cusp 
slope value for the data is a = —1.95 ± 0.02 and normalization 
n = 2.5 ±0.7. 



parameters is doubled. In practice all the data points are 
now inside the one-a Poisson noise range of the model. Pop- 
ulation I parameters are the same as in the single popu- 
lation model. However, the second population parameters 
have much larger relative error range than population HI 
parameters, indicating large uncertainty. This reflects the 
small size of the lens sample. Additionally, the second popu- 
lation cusp slope has such a value that is much more affected 
by the perturbations than the cusp slope value for popula- 
tion HI lenses. 



5.3 Is the second population real? 

Regardless of the Poisson noise limits, we are very concerned 
if the observed feature at o?csl ~ —1.5 is nonetheless noise 
induced. The statistical fit with such a small sample com- 
pared to the number of fitted parameters should be regarded 
with high suspicion. Additionally, the signature of the sec- 
ond population is relatively weak when considering the error 
limits for the fitted parameters of the second population. 

Therefore, we set up a Monte-Carlo experiment in or- 
der to find out the probability that the second population 
hypothesis could be rejected. We generated 100 lens cata- 
logues with 40 lens systems by randomly sampling a mock 
catalogue of 4000 lens systems. The mock catalogue was cre- 
ated using a single lens population with a cusp slope value of 
a = —1.95 and accounting for the lens ellipticities, the sub- 
structure and the source variability. These generated cata- 



Figure 13. A dual population fit of a re-normalized theoreti- 
cal probability function for the CSL value distribution from an 
observational lens catalogue. Jagged histogram and correspond- 
ing diamond symbols with error bars present histogram data for 
ctcsL values calculated from the data. Error bars correspond to 
measurement errors in photometry. Thick solid curve represents 
fitted CSL-distribution, with a (dashed line) and 2<r (shaded re- 
gion) outlining the Poisson shot noise error region. The fitted cusp 
slope value for the data is a H i = —1.95 ± 0.02 and normalization 
n H1 = 2.5 ± 0.7 for the first population. The second population 
has corresponding determined values a H 2 = —1.49 ± 0.09 and 
n H2 = 0.5 ±0.3. 



logues were fitted with a single population model. The prob- 
ability to reject the second population hypothesis is equal 
to the probability to get as strong signature of the second 
population as in the observed lens catalogue. 

In the observations, the second population feature ex- 
ceeds the 2cr noise limit by factor 1.8 at ocsl range 
[—1.6, —1.4]. In order to interpret the noise induced feature 
as the second population, we set a condition that the data- 
set must contain a data point that exceeds the 2a level by 
similar factor at the same range. We found one such cata- 
logue out of one hundred, thus we estimate that the proba- 
bility for rejection of the second population hypothesis is at 
order of one per cent. 

On the other hand, the observational lens data contains 
two data points above the 2a noise level. If we count lens 
samples that have two data points exceeding the 2a thresh- 
old in o?csl range [—1.6, —1.4], we get a rejection probability 
of two per cent. If we add the factor of 1.8 excess condition, 
the probability goes down below one per cent. 

Hence, we estimate that the rejection probability for the 
dual population hypothesis is roughly ~ 1 per cent. 



6 DISCUSSION AND CONCLUSIONS 

Our analysis is based on a geometric measure of a lens sys- 
tem - how much cuspiness does an axially symmetric lens 
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need in order to produce observed magnification ratio, i.e. 
the CSL value. In a sense, analogously to the curvature char- 
acterizing properties of space, it can be thought as a gauge 
characterizing localized properties of a lensing event. There 
is no need to define the exact global properties of the lensing 
potential in order to have a meaningful measure. 

The CSL measure can be derived for axially symmetric 
lenses, and as a result, all the coefficients containing infor- 
mation on cosmology, cosmological angular separations of 
the lens and the source, lens mass and concentration are 
cancelled out. Only preserved quantities are the cusp slope 
and the image magnification ratio. 

If we assume that the observed lenses have completely 
random orientations, thus neglecting effects from the mag- 
nification bias, a theoretical distribution of the CSL values 
can be calculated for a population of lenses. In that case, 
the properties of the lensing profile(s) determine the shape 
of the distribution. When we calculate the CSL value dis- 
tributions from the real life data, we make an additional 
assumption that the observed image fluxes are determined 
only by the lens magnification of a constant flux source. 

This chain of assumptions has its weak points, and we 
have tried to address most of them in section [4] Here we 
are going to do a brief summary of these, and discuss few 
unaddressed ones. 

The first unaddressed issue is at the derivation of the 
magnification ratio equation (|13|l . The equation is based on 
the inner part of the lens equation that does not behave cor- 
rectly when the radial image coordinate at the lens plane ex- 
ceeds the value k B ■ Usually, this happens outside the strong 
lensing region, but if a < —2.0 and the mass of the lens or 
the concentration is high enough, k B can be located inside 
the strong lensing area (/cb *\ k cr 2*). 

This asymptotic behaviour of the lens equation was ac- 
counted for when the mock catalogues were created with the 
Monte Carlo method. The method works correctly with the 
mock data because the deviation from the correct value at 
k > k B happens very slowly with the increasing radial image 
coordinate. Usually, difference becomes non-negligible only 
for images far into the weak lensing region even if k B were 
inside the strong lensing area. Significant deviation requires 
unrealistically high mass or concentration for the halo. 

The next issue is the error factors considered in section|4] 
- lens ellipticity, lens potential substructure, and source vari- 
ability. Our earlier treatment was focused on examining the 
perturbations in the recovered cusp slope value. The anal- 
ysis showed that the most severe disturbations originated 
from the lens ellipticities, although those were at manage- 
able level when a < —1.4 with random ellipticities. 

The aforementioned factors have little or no effect be- 
cause the analysis is restricted to the double image lenses by 
definition. At the double image region, deviations from the 
axially symmetric model are at their weakest. Differences 
become severe nearby the critical curves (at the origin I ~ 
and at the border of the strong lensing region I ~ imax), 
where the magnification diverges. Images nearby the critical 
curves produce either the Einstein rings (I ~ 0), quad lenses, 
or triple image lenses (I < Z max ). Hence, the real life lensed 
images nearby critical curves are easily identified, and they 
can be excluded from the observational sample. 

The similar reasoning can be used to counter arguments 
for effects from the sheared potential, micro-lensing, and the 



magnification bias mentioned earlier. The bias cannot be 
very strongly present in the sample because the lens systems 
with excessive magnifications are excluded from the analysis. 

Because the CSL analysis method is based on the flux 
ratios of the images, the extinction effects from intergalac- 
tic matter that affect both images in a similar way, can not 
change the results. In principle, considerable amount of dif- 
ferential extinction present at the lens object can corrupt the 
CSL value. However, the lensed images are very often well 
outside the visible lens object, which makes it improbable. 

At this point we can conclude that the same reason 
that makes double image lenses unattractive targets for the 
substructure or the micro-lensing studies gives them ideal 
properties as a probe for the macro-lens profile. They are 
resilient against perturbations. 

The main results from our study are the two halo profile 
populations. We get cusp slope value q H i = —1.95 ± 0.02 
and normalization n H i = 2.5 ± 0.7 for population HI and 
for population H2 correspondingly a H 2 = —1.49 ± 0.09, 
ji H 2 = 0.5 ± 0.3. The normalization values indicate that 83 
per cent of the lens objects belong to population HI and 
17 per cent to population H2 with 4 per cent error margins 
in both values. This suggests that, on average, one sixth of 
the lens objects in our sample belongs to the second popula- 
tion with shallower cusp slope. We performed a Monte Carlo 
study in order to find out the rejection probability for the 
second population hypothesis. The result suggests that the 
probability for the second population signature being noise 
induced is at order of ~ 1 per cent. 

As our Monte Carlo testing in section|4]pointed out, the 
profiles with shallower cusp slopes are more easily perturbed 
by error sources. These deviations have tendency to make 
the acquired cusp slope steeper than it really is. Although 
this is partly reflected by rather wide error range in the 
population H2 values, our referred population H2 results 
should be interpreted with caution. The cusp slope value 
q H 2 can be alleged as an estimate for the lower limit of the 
population H2 profiles. 

The foremost alternative explanation for the H2 popula- 
tion is that they are actually unresolved triple image lenses. 
However, this scenario has several weaknesses: the fused pair 
of images should possess a strong autocorrelation in bright- 
ness variations that is equal to the time delay between the 
images. These unresolved lenses would increase the total 
number of the triple image lenses and their corresponding 
statistical frequency to such level that would contradict the 
acquired cusp slope for the main population HI. If the H2 
lenses were unresolved triple image lenses in HI, the share 
of triple image lenses would be 11 out of 100, and the main 
population should have profile a « —1.89 according to the 
equation (1241) . This is well off from the determined profile 
value for the main population. Furthermore, although the 
H2 population of halos is smaller, they produce distribu- 
tion that is consistent with the model. This should not be 
the case if they were triple image lenses - their distribution 
should be more random. 

According to the conventional wisdom acquired from 
the cosmological N-body simulations, the universal cusp 
slope for the dark matter haloes is around a ~ —1.0... — 1.5, 
depending on the author. However, observations on rota- 
tion curves of normal and LSB galaxies suggest even higher 
values, exceeding a = —1.0. Our results contradict most of 
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these results. The lens den sity profiles has also been stud- 
ied bv lGavazzi et all l|2007l ). They used a power law profile 
with joint weak and strong lensing modelling, and confirmed 
the SIS lens population (HI), but did not observe the shal- 
lower profile population H2. The weak lensing does not have 
very good resolution, and the pure power law profiles do 
not account the projected halo component outside the scal- 
ing radius r s realistically. Additionally, they study average 
properties of 22 lens halos. According to our results, three 
or four lenses in their sample should have shallower profiles. 
The second population signature can easily be lost into the 
residual noise when considering average properties of their 
lensing sample. 



In tere stingly, Gustafsson. Fairbairn fe Sommer-Larsenl 
(|2006h andU omano-Di'az et al. have performed cos- 

mological N-body simulations with baryonic matter com- 
ponent. They have found similar populations for the halo 
profiles as we do in our lens study. In the first population, 
the halo profile is made steeper by the rich baryonic con- 
tent, down to the value a ~ —1.9, which is close to the 
isothermal value. They also discovered a second population 
of haloes that were poor in baryonic content, with cusp value 
q ~ —1.3. It would seem, that although baryonic matter 
density is much lower than density of dark matter, it can still 
have non-negligible effects on the evolution of dark matter. 

Our numerical values are also supported by the addi- 
tional results discussed in appendix and the section 12.11 
Flat cores a > — 1 are very poor in strong lensing, and shal- 
low profiles a > —1.5 should produce excessive number of 
triple image lenses. 

The CSL concept has clear advantages. The method 
uses more realistic lensing model inspired by N-body simu- 
lations. The only required data is sampled double image lens 
photometry. The method can separate intrinsic properties of 
the lensing halos from the extrinsic lensing conditions. Thus 
no information on the lens and the source redshifts, the lens 
masses or concentrations, cosmology, or theoretical distribu- 
tion functions characterizing the lens and source populations 
is needed. This removes many uncertainties that are usually 
present in the statistical lensing studies. 

The method has also disadvantages. When general 
properties of dark matter haloes are considered, there is a 
strong selection bias in our study. Our results concern only 
those haloes that are able to do strong lensing, i.e. haloes 
with cusp slopes steeper than a ~ — 1. Generally, the lens 
objects are massive elliptical or early type galaxies residing 
at the core of a small galaxy cluster which would incorporate 
the dark matter halo of the whole cluster into the lensing 
potential. These kinds of objects are very different from the 
late type spiral galaxies or LSB-galaxies, which usually can- 
not perform strong lensing at all. Additionally, there are 
relatively few observed double image lenses at the moment, 
and the CSL method requires a hand-picked sample of rela- 
tively unperturbed lens systems exhibiting axial symmetry. 

New lens systems are found all the time. Certain pro- 
portion of these are double image lenses. Larger data-set 
produces better results, making the CSL analysis more ac- 
curate in the future. 
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APPENDIX A: CUSP SLOPE LIMIT FOR 
STRONG LENSING 



Our previously 
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publis hed lensing 
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theory in 
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density fcproduced by the GNFW mass profile with 
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where constants a and b depend only on the cusp slope value 
a of the profile. The radial coordinate x on the lens plane 
and similar coordinate y on the sour ce plane are defined 
as in ISchneidc r. Ehlers fc Falcol fl992). Corresponding un- 
normalized lens equation is (a/-l) 
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Here po = Apc,r a /T,„ is the normalized projected density. 
This lens equation has its Einstein radius at y = 0, i.e. 
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Here the scale radius is 
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where M is the total mass of the lens and ci the profile 
concentration. The central density is 
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where 



In the previous equations, the critical density p CI is defined 

as 
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and the critical surface density is 
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In equations (|A7|1 and ()A8|) the Hubble constant is Ho, the 
gravitational constant is G and the speed of light is c. Sub- 
scripts L and S refer to the lens and the source, and sym- 
bol D to corresponding angular cosmological distance and 
2 to redshift. The total energy density of the universe is 
parametrized with Q m , Or and SIa (matter, curvature and 
cosmological constant). 

At the limit where the lens cannot produce multiple 
images, the Einstein radius is reduced to zero and becomes 
imaginary. Thus, a condition for strong lensing can be ac- 
quired from the equation (|A3|I as 
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This equation can be solved numerically with the definitions 
presented above, but it can also be approximated with 



b^ (a+iy 1 {a + 2Y 



and 
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with good accuracy when a > — 1. When these approxima- 
tions are inserted into the condition (|A9[) . we can solve the 
maximum cusp slope for the strong lensing as 
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